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List of symbols and notations 



Generally, upper case letters will be used to represent matrices, 
lower case letters to represent elements of matrices. However, for 
the sake of brevity, it is sometimes preferable to refer to matrices 
while speaking of their elements. For instance, instead of saying: 
"the elements of matrix (I-FT) have poles on the left half plane", 
the words "the elements of matrix" can be suppressed without causing 
ambiguity. Which elements of the said matrix are involved should be 
clear to the reader who follows the reasoning in a continuous fashion. 
The following symbols will be used in the sense as indicated: 

P: plant transfer function matrix 

G: open loop compensating filter transfer function matrix 

C: cascade controller transfer function matrix 

F: overall feedback controller transfer function matrix 

H: inner feedback controller transfer function matrix 

U: system input column matrix 

X: plant input column matrix 

Y: system output column matrix 

E: error column matrix 

T: overall transfer function matrix 

I = identity matrix 

s = Laplace operator 

a, b, a, f$ ... = coefficients of polynomials in s 

m = number of outputs of multivariable plant 

n = number of inputs to multivariable plant 
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But m, n are also used as stated below 



m, n, 

LHP 

RHP 



v, w = orders of polynomials in s 
= left half of the s plane 
= right half of the s plane 
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ON THE REALIZATION OF LINEAR MULTIVARIABLE CONTROL SYSTEMS 

Chapter One 

Multivariable control systems 



I-ls Introduction: 

The problem of synthesis of a linear multivariable control 
system consists in selecting a controller or controllers and a 
compensation scheme such that a set of specifications on transient 
and steady-state responses, disturbances and plant parameter variations 
be satisfied. Furthermore, the controller or controllers selected must 
be realizable in the physical sense. 

A number of papers in the past have treated the synthesis problem 
and are mostly concerned with the satisfaction of specifications. 

Little attention has been paid to conditions concerning physical 
construction of the controllers. Recently R.J. Kavanagh* gave a 
detailed account on his study of realizability of non-interacting 
control starting from a given plants but his problem was "mathematical" 
realizability, dealing with the existence or non-existence of a 
solution, and not with the physical construction of the controllers. 

The purpose of this report is to present a method for recognizing 
whether or not, with a given plant and a given specification, the 
controller or controllers necessary to compensate the system can be 

i * 

built from physical components. It is also shown how a choice is 
possible among different compensation schemes to suit the problems of 
availability of components, reliability of components, cost, ruggedness, 
simplicity, or merely designer's taste. 
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1-2: Matrix representation of multivariable control systems: 



A multivariable control system, just like a single variable 
system, is made up of a plant, or actuator, and one or several 
controllers, or compensators, interconnected in a certain fashion. 

Each of these units will have several inputs and several out- 
puts. It is then convenient to represent them with matrices, each 
element of which will be a transfer function relating an input to 
an output. The (i,j)th element relates the j£h input to the ith 
output. Figure 1 represents a system with plant P and one controller C 
connected in cascade. U, X, Y are column matrices. U represents 
the system inputs; X the plant inputs, Y the system outputs, which 
are the same as the plant outputs. 

Let n be the number of plant inputs and m the number of plant 
outputs, ie: 



— — 




— — 


X 1 




y l 


X 2 




y 2 


O 


and Y = 


0 

0 


X 




0 

y 


n 




''m 


— _ 




— — 



In practice, we always have m< n. It would make no sense to try to 
control more output variables than the number of input variables to 
which we have access. If m = n, P is a square matrix. If m-^ n, for 
convenience in the manipulation of matrices, we can add (n-m) arti- 
ficial outputs, which we will monitor and feed back to the system 
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input. By so doing, we add extra rows to the plant matrix P, 
these extra rows having unity elements on the diagonal and zero 
elements elsewhere. As an example, a 3-output, 5-input plant 
matrix would become this square matrix; 



'll 


P 12 


P 13 


P 14 


p 15 


'21 


CM 

CM 

Q* 


P 23 


CM 

a* 


P 25 


'31 


p 32 


P 33 


p 34 


P 35 


0 


0 


0 


1 


0 


0 


0 


0 


0 


1 



Figure 2 shows what modification has been brought to our multi- 
variable system. One can see that the number of system output is 
still 3. The number of system inputs is also 3 since no more than 
3 different inputs are needed to command 3 system outputs. 

As a result of this artificial manipulation, plant matrices 
and controller matrices can always be made square. Then only square 
matrices shall be dealt with. 

1-3; Interaction and non-interactions 

There have been several attempts to define what is meant by 

12 3 

non-interaction in a multivariable system * ’ . The simplest but 

also most logical one seems to be the following; non-interaction 
is said to occur when there is a one-to-one relationship between 
system inputs and system outputs. In other words each system input 
controls one definite output and that one only. Conversely, when- 
ever one input not only controls its own output but also changes 
one or more other outputs, the system is said to be interacting. 
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Ideally, non-interaction is the goal in the design of multi- 
variable systems. But from an engineering viewpoint, this is not 
always physically realizable, nor is it always desirable, for a 
number of reasons, practical, economical, commercial, among others. 

It is then reasonable to proceed as follows: attempt to design 

the system for non-interaction, and build it if it is judged satis- 
factory for the various reasons mentioned above. If non-interaction 

f-.CCK • ( 

is unrealizable or impractical, then go to the next-to-best solution, 
which is a system with some minimum amount of interaction, or with 
a specified amount of interaction as the case may be, until a 
compromise is reached between satisfactory performance and physical 
and economical realization. 
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Chapter Two 



Compensation design 

II- 1; Selection of a compensation schemes 

In the current literature, a current practice has been to 
arbitrarily select a compensation scheme and solve the problem for 

that particular scheme . Some writers assume one cascade and one 

/| 

feedback controller , others restrict their treatment to a system 
with only a cascade controller^’ Another assume® a feedback 

t 

controller only, and goes as far as restricting himself to a 
diagonal plant matrix^. 

It is true that some compensation schemes are more commonly 
used than others, but it is also true that some other scheme may 
serve a particular purpose better, or suit existing conditions 
better in a given situation. Consequently it is preferable for 
the designer to gain an insight to several compensation schemes 
and try to select the most suitable. 

Figure 3 shows 3 possible structures using both cascade and 
feedback controllers * 

While schemes no 1 and 3 seem to be more general and 
reasonable, in some cases the designer may have to use scheme no 2, 
for example when end-to-end feedback for the overall system is 
impractical or impossible due to geographical conditions 0 

II-2g Filter designs 

In order to avoid restricting the design to any particular 
arrangement, it is proposed that synthesis be based on the simplest 
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system configuration as indicated on Figure 4, where G represents 
a multivariable cascade filter. No feedback of any kind is involved. 

For a given desired overall response, ies a specified overall 
transfer function, say T, we haves 



or 



PG = T 



G 




T 



There is no difficulty in obtaining the expression of G. The 
operations in (1) are matrix operations, hence order is important 
since generally PG 4 GP. 

Note that, as stated before, design was started off by assuming 

( 

non-interaction, therefore T, the overall transfer matrix is diagonal, 
and each diagonal element t^ is the desired transfer function of the 
ith variables 




In the next step, filter G is split into cascade controller C 
and/or feedback controller (s) F and H. Several possibilities then 
occur 

- if none of the C, F, H controllers thus obtained are physically or 
economically realizable, it can be concluded that non-interaction 
is not possible. One must then go to the next best solution, which 
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is the scope of another report to follow 

- if only one of the above schemes is realizable, there is no choice 

but to take it 

- if more than one schemes are realizable, the designer has the 
privilege of choosing one of them in the light of commercial and 
economical conditions, availability of components, reliability, 
ruggedness or his own personal taste. 

The design based on this GP configuration has 2 further advantages; 

- in cases where it is necessary to design the system with inter- 
action, it is easier to see which element is to be redesigned in 
order to make one particular interaction, say the (i, j)th, meet 
the limits specified 

- as will be shown later in this report, examination of G alone is 
sufficient to determine whether or not the compensating schemes 
of Figure 3 are physically realizable. This permits the designer 
to avoid wasting further time on an impossibility. 

In the compensation design using filter G, it is convenient 
to have a set of formulas for going from G to C, F and/or H. 

For scheme no 1 * 

Specifications are given in the form of an overall transfer 
matrix T, or may be converted into such a matrix. The writer 
hopes to undertake a more detailed study of this matter in the 
next paper. By definition of T; 





Y = 


TU 


(3) 


But from figure; 


Y = 


PX 


(4) 




X = 


CE 


(5) 




E = 


U - FY 


(6) 
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Again the order of operations is important since matrix 
quantities are involved „ 



From (4), (5) and (6): Y= PC(U-FY) 



f 


or : [i + PCF] Y = PCTJ 


(7) 


where I is (nxn) identity matrix. 
Premultication of both sides of (7) by Cl + 


pcf]’ 1 yields s 




Y = [l + PCFl^PCU 


(8) 


Comparison of equations (3) and (8) yields % 






T = [i + PCF] _1 PC 


(9) 


But we know that 


T = FG 


(10) 


Equating (9) and 


(10) and premultiplying both sides by [l + PCF] 




PC = [l + pcf] FG 




or 


PC = PG + PGFPG 




Rearranging 


tC [l - FPG] = tG 





Noting that PG = T, 



C = G [i - FT]" 1 



(11) is the relation between G on one hand and C, F on the 
other. This is one equation with two unknowns, C and F. The 
designer is free to choose one, say F, the other variable, C, 
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is then determined. Needless to say, F, that is the elements 
of F, must be chosen so that they are physically realizable. 
It would be helpful if knowledge of G alone can determine 
whether or not a C, F configuration can be realized. This is 
the purpose of the last part of the paper. 

For scheme no 2 - 



T = PG 

Calculations similar to the previous case lead to: 

T = [i + PH] _1 PC 

Equating the 2 equations : PG = [l 4- PH] '‘"PC 

Premultiplying by [i - PH] : PC = PG + PHPG 
4 

Since PG=T, C=G+HT (12) 

Equation (12) relates G with C and H. As before, selection 
of H permits calculation of C or vice versa. 

Note that, in scheme no 1, a possible choice is to take F = I, 
ie: unity feedback without crossfeed. This is a very realistic 
choice since in a control system it is desirable to compare 
each input with its own output on a one-to-one basis. But here 
in scheme no 2, feedback matrix H will generally be non-diagonal 
since there is no one-to-one correspondence between y's and x's. 
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For scheme no 3 - 



T = PG 



Successive reduction of block diagram leads to; 

T = PG = [i + (I + PH)" 1 PCF]" 1 (I + PH)‘ 1 PC 

Premultiplication by [i + (I + PH) _1 PCF] yields; 

PG + (I + PH) _1 PCFPG = (I+PH) _1 PC 
Rearranging; 



PG = (I + PH)" 1 PC Cl - FPG] 
Premultiply by (I+PH); 



(I + PH) PG = PC (I - FPG) 
Postmultiply by (I - FPG)" 1 ; 

(I + PH) PG (I - FPG)" 1 = PC 
or; PC = (PG+ PHPG) (I - FPG)" 1 
or; fC - t (G+HPG) (I - FPG)" 1 



ie 



C = (G + HT) (I - FT) 



(13) 
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Equation (13) relates G with C, F and H. In the most general 
case of unity feedback, F = I and (13) becomes 



C = (G + HT) (I - T)" 1 



(14) 



Selection of one of the 2 unknowns, C and H, determines the 
other. This is the essence of the designer’s task. 

Note that (13) is a more general expression of (11) . Letting 
H = 0 in (13), one obtains (11). 

The following example shows how by use of equations (11) (12) (13), 
a sensible choice between the above compensating schemes is possible. 
Example - Given the plant to be controlled 
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1 A 
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s+20 
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1 


^s+20 


s+3 J 



It is desired that the transient responses y^, y^ have no 
overshoot and have a short rise time as expressed by the 

4 

overall transfer functions t,. = t„ 0 = — . 

11 22 (M-2) 2 

First obtain the transfer matrix for filter Gs 



G = FT 
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